Accurate calculations of Sr properties for a high-accuracy optical clock 
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We have carried out calculations towards the goal of reducing the inaccuracy of the Sr optical 
atomic clock to lxlO~^^ and below. We calculated a.c. polarizabilities of the 5s^ ^So and 5s5p ''Pq 
clock states that are important for reducing the uncertainty of blackbody radiation-induced fre- 
quency shifts for the ^5*0 — "^Pq clock transition. We determined four low-lying even-parity states 
whose total contribution to the static polarizability of the Pq clock state is at the level of 90%. We 
show that if the contribution of these states is experimentally known with 0.1% accuracy, the same 
accuracy can be achieved for the total polarizability of the Pq state. The corresponding uncertainty 
for the blackbody shift at a fixed room temperature will be below IxlO"'^^. The calculations are 
confirmed by a number of experimental measurements on various Sr properties. 

PACS numbers: 31.15.ac, 31.15.am, 31.15.ap, 32.70.Cs 
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I. INTRODUCTION 

The use of optical lattice-confined neutral atoms for 
the goal of achieving a new level of time-keeping precision 
and accuracy has become widely practiced (see, e.g., [l|, 
0,Slj,Sfi[3,SQ)- In this scheme, ultracold atoms are 
confined in an optical lattice to eliminate motion-related 
systematic frequency shifts. The lattice laser wavelength 
is selected in such a way that the perturbation to the 
clock frequency arising from a.c. Stark shifts induce d by 
the lattice laser for both clock states exactly cancel [10l |. 

One indicator of clock performance is provided by the 
Allan deviation characterizing the fractional frequency 
instability a. For a signal-to-noise ratio given by the fun- 
damental quantum projection noise, the instability can 
be written as: 



a{T) 



1 ^/T 



(1) 



Here Q is the resonance quality factor defined as Q = 
V I ts.v, where v is the transition frequency and IS.v is the 
lincwidth. N is the total number of particles measured 
in a coherent time period of T, and r is the total averag- 
ing time. According to Eq. ([1]), atoms with the highest 
quality factors Q are preferred for a new generation of 
time and frequency standards. 

The highest Q's are currently obtained for a narrow 
transition in the optical domain [ll|. In particular, the 
highly forbidden ^5*0 -^ ^Pq transitions in divalent atoms 
offer excellent possibilities for attaining a new level of 
precision and accuracy for time-keeping. One of the 
most promising candidates is ^^Sr for which Ai/ « 1.2 
mHz [H [H til, yielding a potential Q > 10^^. In 
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a recent paper [g, a systematic uncertainty evaluation 
for a neutral Sr optical atomic standard was reported 
at the 10^^^ fractional level, surpassing the best evalu- 
ations of Cs fountain primary standards. The dominant 
systematic frequency correction and uncertainty in that 
work arose from the room temperature blackbody radia- 
tion (BBR). The fractional frequency shift, {Svebr/i^oI, 
caused by the BBR is proportional to the differential 
static polarizability of the two clock states. For the 
5s^ ^5*0 — > 5s5p '^Pn transition in Sr, the BBR shift was 
calculated in Ref. [3 to be equal to 55.0(7) xlQ-i^. The 
1% uncertainty for the BBR shift originates mostly from 
insufficient knowledge of the static polarizability of the 
5s5p "^Pg state, with the most accurate calculation pro- 
vided in [l5|. To further improve the Sr accuracy, it 
is clear that better understanding of the Sr properties is 
needed to give a more accurate determination of the BBR 
shift. The purpose of this paper is to outline a clear path 
to achieve this goal. To improve the clock accuracy signif- 
icantly, it is equally important that a well-characterized 
homogeneous BBR environment surrounds the Sr atoms 
in future experiments. 

The improvement of the Sr clock accuracy requires a 
more accurate determination of the differential static po- 
larizability of the 5s^ ^Sq and 5s5p ^Pq states. Note that 
the static polarizability of the ground 5s^ ^5*0 state is 
known at a sufficiently low uncertainty ~ 0.1% [l^- This 
low uncertainty is made possible by a good knowledge 
of the matrix element \{5s^ %||d||5s5p ^Pf)] [3, where 
the intermediate state 5s5p ^P° contributes to the polar- 
izability of the ^Sq state at the dominant level of 97%. 
Consequently, the outstanding challenge is to reduce the 
uncertainty of the polarizability of the ^Pq state to a sim- 
ilar level. Even sophisticated modern relativistic meth- 
ods of atomic calculations cannot provide such accuracy. 
For this reason, a solution to this problem must combine 
theoretical and experimental approaches. We show that 
four specific intermediate states have a combined contri- 
bution to the total static polarizability of the 5s5p '^Pq 



state at the level of 90%. When the contributions from 
these four states are determined from experimental data 
at 0.1% accuracy and the contributions of all other dis- 
crete and continuum states are known at the level of 1-2% 
from calculations, then the final 0.1% uncertainty for the 
polarizability of the bsbp ^P§ state will be achieved. This 
strategy is the focus of this paper. 

The paper is organized as follows. In Section II we 
briefly describe the method of calculations. In particu- 
lar, we discuss the construction of basis sets and solving 
the multiparticle Schrodinger equation. In Section III we 
discuss the blackbody radiation effect and present the re- 
sults of calculations for the low-lying energy levels, a.c. 
polarizabilities, transition rates, and oscillator strengths, 
and we analyze the results obtained. Section IV contains 
concluding remarks. Atomic units (h = \e\ = m = 1) are 
used throughout the paper. 



II. METHOD OF CALCULATIONS 

The most complex problem in precise atomic calcu- 
lations is associated with the necessity to account for 
three types of electron correlations, i. e., valence- valence, 
core-valence, and core-core correlations. The former are 
usually too strong to be treated perturbatively, while 
the other two types of correlations cannot be treated ef- 
fectively with non-perturbative techniques, such as the 
multi-configurational Hartree-Fock metho d fl7| or the 
configuration interaction (CI) method [18l.ll9|. 

Therefore, it is natural to combine the many-body 
perturbation theory (MBPT) with one of the non- 
perturbative methods. In Ref. [201, it was suggested 
to use MBPT in order to construct an effective Hamil- 
tonian for valence electrons. After that, the multiparti- 
cle Schrodinger equation for valence electrons is solved 
within the CI framework. Doing so allows us to find the 
low- lying energy levels. Following the earlier works, we 
refer to this approach as the CI-I-MBPT formalism. 

In order to calculate other atomic observables, one has 
to construct the cor resp onding effective operators for va- 
lence electrons [2l|, [2^ [23[. These operators effectively 
account for the core-valence and core-core correlations. 
In particular, to obtain an effective electric-dipole opera- 
tor, we solve random-phase approximation (RPA) equa- 
tions, thus summing a certain sequence of many-body 
diagrams to all orders of MBPT [Ili [13, [l^. The RPA 
describes shielding of an externally applied field by core 
electrons. Small corrections due to, for instance, normal- 
ization and structural radiation are omitted. 

In the CI-I-MBPT approach, the energies and wave 
functions are determined from the time-independent 
Schrodinger equation 



iJoff(£^„)$„ = En'^r. 



where the effective Hamiltonian is defined as 



H,siE) = Hfc + S(^). 



Here -ffpc is the Hamiltonian in the frozen core approxi- 
mation and S is the energy-dependent correction, which 
takes into account virtual core excitations. The opera- 
tor E completely accounts for the second-order pertur- 
bation theory over residual Coulomb interaction. De- 
termination of the second-order corrections requires cal- 
culation of one- and two-electron diagrams. The one- 
electron diagrams describe an attraction of a valence 
electron by a (self)-induced core polarization. The two- 
electron diagrams are specific for atoms with several va- 
lence electrons. The number of the two-electron dia- 
grams is very large and their calculations are extremely 
time-consuming. In the higher orders the calculation 
of two-electron diagrams becomes practically impossible. 
Hence, it is more promising to account for the high- 
orders of the MBPT indirectly. One of such methods 
was suggested in Ref. [2g|, where it was shown that a 
proper choice of the optimum initial approximation for 
the effective Hamiltonian can substantially improve the 
agreement between calculated and experimental spectra 
of many-electron atom. 

We consider Sr as a two-electron atom with the core 
[ls,...,4p^]. The one-electron basis set for Sr includes Is- 
14s, 2p-14p, 3d-13d, 4/-13/, and 5g-9g orbitals, where 
the core- and 5s-7s, 5p-7p, and M-dd orbitals are Dirac- 
Hartree-Fock (DHF) ones and all the rest are the virtual 
orbitals. The orbitals ls-5s were constructed by solving 
the DHF equations in the V'^ approximation, i.e. the 
core and the 5s orbitals were obtained from the DHF 
equations for a neutral atom (we used the DHF com- 
puter code [23|). The 6s-7s, 5p-7p, and 4d-6d orbitals 
were obtained in the V^~^ approximation. That is, the 
ls-5s orbitals were "frozen" , one electron was transferred 
from the valence 5s shell into one of the orbitals specified 
above, and the corresponding one-electron wave function 
was found by solving the HFD equations. We determined 
virtual orbitals using a recurrent procedure similar to 
Ref. J^] and described in detail in ^, 221 ■ 

Configuration-interaction states were formed using 
these one-particle basis sets. It is worth emphasizing 
that the employed basis set was sufficiently large to ob- 
tain numerically converged CI results. An extended basis 
set, used at the stage of MBPT calculations, included Is- 
21s, 2p~21p, 3d-20d, 4/~17/, and 5g~13g orbitals. 



III. RESULTS AND DISCUSSION 

A. Calculation of energies 

Solving the multiparticle Schrodinger equation we find 
low-lying energy levels and their respective wave func- 
tions. In Table U we present the calculated energies of 
the low-lying states for Sr and compare them with exper- 
imental data. As is seen from the table we focus mainly 
on the energy levels with J = 1. This is due to the 
fact that we are interested in calculation of the electric 
dipole-dominated a.c. polarizabilities for the clock states 



with total angular momentum J — 0. Only intermedi- 
ate states with J — 1 contribute to these polarizabilities. 
The energy level diagram of these states is given in Fig.[T] 
The energy levels were obtained in the framework of the 
conventional configuration-interaction method as well as 
using the formalism of CI combined with the many-body 
perturbation theory. Using the CI method alone, the 
agreement of the calculated and experimental energies is 
at the level of 5-10%. The combination of CI and MBPT 
improves the accuracy by approximately an order of mag- 
nitude. 
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FIG. 1; (Color online) Low-lying energy levels and transition 
wavelengths of atomic Sr, relevant for the determination of 
polarizabilities of the 5s^ ^5*0 and 5s5p ^Pg clock states. 



B. Blackbody radiation 

It is well known that BBR-related clock frequency 
shifts arise from perturbations of atomic energy lev- 
els by weakly oscillating thermal radiation. For the 
-^^o -^ ^P§ clock transition, both atomic states involved 
are perturbed. Thus, the net BBR effect is the differ- 
ence of the individual BBR shifts for the two states, 
<5wBBR = 5Ebbk{^P§) - <5£'bbr(^S'o). The expression for 
the SE-Q-mig) of a 5 state can be given by the following 
formula llSlI 



5E 



BBR 



10 



(2) 



where a ~ 1/137 is the fine structure constant, T is 
the characteristic temperature of the BBR environment, 
aEg (0) is the electric-dipole static polarizability of the g 
state, and 77 represents a "dynamic" fractional correction 
to the total shift. As was shown in Ref. [I^, rj is negligi- 
ble for the ^S'o state but contributes to (J-Ebbr of the ^Pq 
state at the 2.7% level. This is primarily due to the fact 



TABLE I: Energy differences (in cm ^) with respect to the 
ground Ss'^ ^So state for the low-lying energy levels of Sr. 



Config. 


Level 


CI 


CI+MBPT 


Experiment [29] 


5s4d 


'D, 


19571 


18076 


18159 


5s4d 


'D, 


19587 


18141 


18219 


5s4d 


'Ds 


19617 


18254 


18319 


5s4d 


'D, 


20166 


19968 


20150 


5s6s 


^5i 


27488 


29019 


29039 


5s5d 


'D, 


33358 


35060 


35007 


5p^ 


'Pi 


33511 


35326 


35194 


5s7s 


's. 


35695 


37429 


37425 


5s6d 


'D, 


37985 


39725 


39686 


5s5p 


'Po 


12490 


14241 


14317 


5s5p 


■.ipo 


12663 


14448 


14504 


5s5p 


'Pi 


13022 


14825 


14899 


5s5p 


ipo 


20832 


21469 


21698 


5s6p 


3po 


32110 


33814 


33868 


5s6p 


ipo 


32487 


34105 


34098 


4d5p 


'D°, 


36699 


36189 


36264 


4d5p 


■ipo 


36944 


37213 


37303 


5s7p 


ipo 


37275 


38927 


38907 


5s7p 


■ipo 


37939 


39377 


39427 



that the ^Pq state has a transition frequency to a nearby 
BsAd ^Di state in the infrared, as shown in Fig. [1] 

The electric dipole static polarizabilities of the 5s^ ^5'o 
and 5sbp "^Pq states were calculated in Ref. [l^ using 
the same method, but a different basis set, as pre- 
sented in this work. The corresponding relative fre- 
quency shift of the clock transition was determined to 
be Idi^BBBjj^ol = 55.0(7) x 10"^''. The shift uncertainty 
of AISvbbr/i^o] = 0.7 X 10^^^ results directly from the 
1% uncertainty attained in Ref. [1^] for the static po- 
larizability A[a3po(0)] of the 5s5p "^Pq state. The small 
size of rj ensured that no additional uncertainty from the 
dynamical correction contributed at this level. 

An equally important source of uncertainty in the 
actual BBR shift is the knowledge and control of the 
blackbody environment at room temperature, T. From 
Eq. ^, the shift uncertainty A[SEBBR{g)] originating 
from the uncertainty in the BBR environment AT is 

A[SEBBR{g)] - -^(a^)3T3a£,(0)[l-Kr;]Ar, (3) 

At room temperature, measurement of the BBR envi- 
ronment at the uncertainty level of AT = 1 K leads to a 
fractional frequency shift uncertainty of 7.5 x 10"^^ 0. 
The combination of this uncertainty in quadrature with 
that resulting from A[a3po(0)] yields a 1 x 10"^^ total 
BBR uncertainty, which currently limits the accuracy of 
the Sr optical clock. 

To further improve the Sr accuracy, the total BBR 
uncertainty must be reduced. This requires solving two 
main problems: i) measuring and controlling the black- 
body environment to much better than AT = 1 K; ii) 
determining the differential static polarizability to bet- 
ter than 1% uncertainty. The first of these requires ad- 



ditional care and design in the experimental apparatus. 
While the temporal stability of the BBR temperatm'e is 
more or less straightforward to achieve, the experimental 
difficulty originates from achieving spatial homogeneity 
of the BBR temperature over various functional areas of 
the vacuum chamber housing the atomic sample. Typ- 
ically, a large fraction of the Air sterradians of solid an- 
gle around the atoms consists of glass viewports to ac- 
commodate the optical access needed for various atomic 
manipulations (laser cooling and trapping, loading into 
an optical lattice, state preparation, etc. ). These areas 
are more difficult to precisely temperature stabilize than 
the remainder of the solid angle typically composed of 
metallic vacuum chamber. Experimentally, we observe 
that different parts of the Sr vacuum apparatus at JILA 
can vary by as much as 1 K. Furthermore, at the high- 
est clock accuracy level, it is important to account for 
the effect of the transmissivity of glass viewports for vis- 
ible and infrared radiation from the ambient room on the 
blackbody environment seen by the atoms. One approach 
to reducing the uncertainty of the BBR shift is to sur- 
round the atomic sample in a cryogenically-cooled shield 
[301 . Doing so reduces both the magnitude and thus the 
uncertainty of the BBR shift. Another approach is to 
enclose the atoms in a chamber closely resembling the 
blackbody cavities used for thermal radiation metrolog- 
ical standards |31|. For example, the optically-confined 
atoms can be transported in a moving lattice from a 
main chamber to a smaller, blackbody cavity [a|. By 
careful temperature control of this small cavity made of 
highly-thermally-conductive material, excellent temper- 
ature homogeneity can be maintained. The very limited 
optical access (for lattice laser and clock probes) enables 
the effective emissivity of the cavity interior to be very 
close to unity. To reach the 10~^^ clock uncertainty, the 
BBR environment must be known at the part per thou- 
sand level at room temperature, corresponding to a BBR 
temperature accuracy at the 100 mK level. 

The differential polarizability must also be carefully 
determined to higher accuracy. In the case of cesium, 
a well-controlled d.c. electric field has been used to in- 
duce a clock shift and determine the differential static 
polarizability [33] ■ As well, some atom interferometric 
techniques may hold promise for directly measuring the 
differential polarizability at better than the 1% level 33]. 
Here we address the improved determination of the differ- 
ential polarizability based on atomic structure measure- 
ments. The uncertainty of the differential static polariz- 
ability is determined by the uncertainties in the polariz- 
abilities of the two clock states. The static pola rizability 
of the ground state ai5„(0) = 197.2(2) a.u. fTS] is known 
with 0.1% accuracy. Consequently, the task at hand is to 
determine the static polarizability of the 5s5p ^Pq state 
with a similar level of accuracy. This is a key step to- 
wards Sr lattice clock operation at the 10^^^ uncertainty 
level. We now discuss this problem in detail and present 
a possible solution in the following sections. 



C. Calculation of electric dipole a.c. 
polarizabilities of the 5s^ ^So and 5s5p ^Pq states 

Using the wave functions of the low-lying states ob- 
tained as a result of solving the multiparticle Schrodinger 
equation, we are able to calculate a.c. polarizabilities of 
the 5s^ ^So and the 5s5p ^Pj states. As one check of 
the quality of our calculations, we can find the magic 
wavelengths: Aq at which ai5j,(Ao) ~ aspo (Ao) and Ai 
at which Q;i5g(Ai) — aapo (Ai) and compare these val- 
ues against the experimental results. In recent works 
these magic wavelengths were determined with high pre- 
cision to be Ao(% - ^Pg) = 813.42735(40) nm g and 
Ai(% - ^P1{m,j = ±1)) = 914(1) nm for linear polar- 
ization [3J|. Furthermore, our calculation can also be 
checked against a recent measurement of the a.c. Stark 
shift associated with the probe of the (^5*0 ~> ^Pq) dock 
transition itself p\. 

We start with a brief description of the method used 
to calculate the electric dipole polarizabilities. The equa- 
tion for the a.c. electric dipole polarizability of the g state 
can be written in the following form, 

as, [uj) ^ 2 2_^ 



E 



l(gMo|fc)P 

Ek - {Eg + u) 



E 



l(gMo|fc)P 

Ek - [Eg - u) 



2{«s.+c.(0) + aB,-c.(0)} 



(4) 



The two terms in the bottom-line of Eq. Q can be viewed 
as the static polarizabilities of the g state calculated for 
the shifted energy levels oi Eg + u and Eg — uj, respec- 
tively. Thus, our task is reduced to computation of these 
two static polarizabilities. 

Following Refs. [33, [3y| we decompose an a.c. polariz- 
ability into two parts. 



a{u!) = a^^Lu) + a'^i^Lu). 



(5) 



The first term describes excitations of the valence elec- 
trons. The second term characterizes excitations of core 
electrons and includes a small counter term related to 
excitations of core electrons to occupied valence state. 
The core polarizability a'^ was calculated at uj=0 to be 
a'^(O) = 5.4 a.u. [1^. Since a'^ contributes to the to- 
tal polarizability only at the level of a few percent and 
its dependence on frequency is very weak, the value of 
5.4 a.u. can also be used for calculations of the total 
a.c. polarizabilities. This approximation of a constant 
core polarizability over the relevant frequency range in- 
troduces an additional uncertainty of <0.1% to the total 
^Pq polarizability. 

It is worth mentioning that the core is the same for 
the 5s^ ^Sq and the 5s5p ^Pj states. For this reason 
a'ig [uj) « Qfapo {ijj) and we arrive at the following ex- 
pression. 



ns„(w) - azpo{uj) « a1s^{uj) 



*3po 



(c). 



(6) 



TABLE II: Calculated polarizabilities at a few selected optical 
wavelengths. Wavelengths A are in nm, the frequencies uj 
are in a.u. and the electric dipole a.c. polarizabilities of the 
Ss'^ ^5o and the 5s5p ^-Pqi states are in a.u. The polarizability 
of the P° state is calculated for the projection \mj\ = 1 and 
linearly polarized light. 



A 


LU 


a(5s2 %) 


a(5s5p ^PS) 


a(5s5p ^P°) 




0.0000 


197.2 


457.0 


498.8 


698.4 


0.0652 


351.8 


909.2 




805.0 


0.0566 


288.9 


289.3 




813.4 


0.0560 


286.0 


280.5 




902.2 


0.0505 


263.5 




263.4 


914.0 


0.0499 


261.2 




256.1 



The method of calculation for the dynamic valence po- 
larizabilities a"(u;) is described elsewhere (see, e.g., [3a, 
|37|1. Here we only briefly recapitulate its main fea- 
tures. These polarizabilities are computed with the 
Sternheimer [S^ or Dalgarno-Lewis 39] method imple- 
mented in the CI-I-MBPT framework. (Here we denote 
E and RPA corrections as the many-body perturbation 
theory (MBPT) corrections.) Given the g state wave 
function and energy Eg, we find intermediate-state wave 
functions (5'0± from an inhomogeneous equation, 



\Hi 



1 



H,s - {Eg ± uj) 



Y,\k){k\do\9) 



1 



■415)- 



TJeff - {Eg ± uj\ 
Using Eq. (|4]) and 5'4)± introduced above, we obtain 



(7) 



(8) 



where superscript v emphasizes that only excitations of 
the valence electrons are included in the intermediate- 
state wave functions 5^± due to the presence of -ffcff ■ 

In Table |TT] we present the values of the static polariz- 
abilities and the a.c. polarizabilities of the 5s^ ^5*0 and 
the 5s5p ^Pqi states computed for different values of A 
using the CI-I-MBPT approach. As a first step we solved 
an inhomogeneous equation and found the valence parts 
of the polarizabilities. Then the values of the polarizabil- 
ities of the ground state were corrected as follows. We 
used the fact that the intermediate 5s5p ^Pf state con- 
tributes to this polarizability at the level of 97%. Know- 
ing the experimental energy difference {Eipo — Eis„ ) and 
the matrix element \{5s^ ^So\\d\\5s5p ^P^)\ = 5.249(2) 
a.u. extracted from the precise measurement of the life- 
time of the 5s5p ^Pi state [lal, we replaced the theo- 
retical contribution of the 5s5p ^P^ state to the ground- 
state polarizability by the experimental value. Finally, 
we added a'^ term to the valence parts, arriving at the 
values listed in Table HTl 

Starting from the 0.05% uncertainty of the 
|(5s^ ^5011^11535^ ^Pf) I matrix element, we estimated 



the uncertainty of the a.c. polarizability of the ground 
state at the level of 0.1%. In particular, for the static 
polarizability, we obtained 015,-, (0) = 197.2 a.u., in 
perfect agreement with the result obtained in Ref. [1^1 
using a different basis set. 

Experiments [8|, |3J] have determined the magic wave- 
lengths for the iS'o - ^P§ and % - ^P^{mj = ±1) transi- 
tions to be 813.4 nm and 914 nm, respectively. As is seen 
from Table |IT1 the calculations carried out in this work 
give the values of 805 nm and 902 nm for these magic 
wavelengths, respectively. Thus, the agreement between 
theoretical and experimental results is at the level of 1%. 
The behavior of the a.c. polarizabilities of the ^Sq and 
the ^Pq states in the wavelength range from 650 nm to 
950 nm is illustrated in Fig.[2l A large peak at 679 nm for 
the ^Pq state arises from the contribution of the 5s6s ^Si 
state, while a small peak in the vicinity of 690 nm for 
the ^^o a.c. polarizability is due to the contribution of 
the 5s5p '^Pi state. Experimentally, the differential a.c. 
polarizabilities in the form of the clock frequency shift 
induced by the clock probe laser itself is known. With 
a probe laser intensity of 20 mW/cm^, the fractional 
frequency shift was measured -1.5(0.4)x 10~^^ [S]. As- 
suming the same probe laser intensity and the values of 
the polarizabilities for the ^5*0 and ^Pf states obtained at 
698.4 nm (see Table lll| . the calculated fractional shift is 
-1.2 X 10~^^, in a good agreement with experiment. 
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FIG. 2: (Color online) Electric dipole a.c. polarizabilities 
for 5s^ ^5'o (solid line) and 5s5p "^Po (dashed line) states of 
Sr. The polarizabilities are shown as a function of optical 
wavelength A. 

It is worth noting that knowing the precise exper- 
imental values of the magic wavelengths and using 
the fact that aapo (813.4 nm) = aig^, (813.4 nm) and 
«3po.|m^l=i(914 nm) = ai5^(914 nm), we can refine our 
calculation and predict with high accuracy the values of 
the a.c. polarizabilities of the ^Pq and ^Pf states at these 
wavelengths. We obtain aapo (813.4 nm) = 286.0(3) a.u. 
and aapo |„j^|^]^(914 nm) — 261.2(3) a.u., matching the 
polarizabilities of the ^Sq state at these two wavelengths. 



TABLE III: Individual contributions (I.C.) from six low-lying intermediate states to the valence parts of the static polarizability 
a3po(0) = 451.5 a.u. and to the a.c. polarizabilities (in a.u.) at the wavelengths A — 813.4 nm and A — 698.4 nm. cvapo (813.4 

nm) = 280.6 a.u. and aapo (698.4 nm) = 903.8 a.u.. D = l(5s5p Pq ||d||n)| is the reduced matrix element of the electric dipole 



operator d. The row "Total" gives the sum of the contributions for each column. 



D (a.u. 



I.C. to a^pg (0) 
u. " % 



a.u 



I.C. to ajp„ (813.4 rmi) 

% 



I.C. to aSp„ (698.4 nm) 
a.u. ° % 



5s4d Tk 
5s6s ^5i 
5s5d ^A 
5p" ^'Pi 
5s7s ^51 
5s6d ^A 



2.74 


286.1 


63.4 


-31.0 


-11.3 


-22.2 


-2.5 


1.96 


38.3 


8.5 


126.4 


45.9 


705.8 


78.1 


2.50 


44.2 


9.8 


68.3 


24.8 


84.8 


9.4 


2.56 


45.3 


10.0 


68.7 


25.0 


84.1 


9.3 


0.52 


1.7 


0.4 


2.4 


0.9 


2.8 


0.3 


1.13 


7.4 


1.6 


9.6 


3.5 


10.8 


1.2 



Total 



423.4 



93.8 



245.8 



89.3 



868.4 



96.1 



In Table IIIII we present the values of the individual 
contributions from six low-lying even-parity intermedi- 
ate states to the valence parts of the static polarizability 
and the a.c. polarizabilities at A = 813.4 nm and A — 
698.4 nm for the ^Pg state. As is seen from Table |llTl 
for the ^Pq static polarizability, four even-parity states 
[hsM^Di, bs&s^Si, 5s5d^Di, and 5p^ ^Pi) contribute 
at the level of 92%. If the sum of these four contributions 
are determined experimentally with an accuracy ^0.1% 
and the contributions of all the rest of the discrete and 
continuum states are known at the level of 1-2% from cal- 
culations, the overall polarizability of the 5s5p ^Pq state 
can be determined with the accuracy ~0.1%. 

Experimental determination of the four dominant con- 
tributions may be accomplished directly from lifetime or 
transition rate measurements. However, the lifetime data 
should be accompanied by high accuracy branching ra- 
tio measurements. Alternatively, we could measure single 
channel decay directly to the 5s5p ^Pq state. In addition, 
it is possible to constrain the four contributions using 
other spectroscopic data such as the magic wavelength 
Ao and the light shift at 698 nm, which will naturally 
be measured and confirmed by a number of future clock 
experiments. At both wavelength regions the same four 
states dominate the polarizability as in the static case. 
The general strategy would be to use spectroscopic data 
to constrain the most dominant contribution for that spe- 
cific case. 

Since the 5sAd '^Di state dominates the four critical 
contributions to the static polarizability, and has less con- 
tribution to the dynamic polarizability at wavelengths 
of interest (see Table HlH) . it would thus be maximally 
beneficial to measure this contribution directly with an 
oscillator strength measurement between 5s5p ^Pq and 
5s4(i ^Di. Doing so avoids upscaling of the uncertainty 
via error propagation. The high accuracy experimental 
measurement for the 5s5d ^Di state 50] can be further 
improved by monitoring its decay to 5s5p ^Pq directly. 
Then, combining the measurements of Aq and the light 
shift at 698 nm would allow us to constrain contributions 
from both 5s6s ^Si and 5p^ ^Pi , permitting constraining 
a3po(0) at the level of ais„(0). 



The magic wavelength, Ai, could in principle also aid 
the constraint. However, additional sensitivities of the 
bs^p^Pi polarizability from J = and J = 2 even-parity 
states are large and complex. Furthermore, the vector 
nature of this state requires careful experimental control 
of light polarization. 

In the following section we present the results of the 
theoretical calculation of transition rates and oscillator 
strengths most relevant to the ^Pq and -^5o polarizabili- 
ties and compare them to existing experimental and the- 
oretical data in the literature. 



D. Transition rates and oscillator strengths. 

The transition rate {W) and the oscillator strength (/) 
for a transition from an initial state \^' J' L' S') to a final 
state (7 JLS"! can be represented as (see, e.g., |40| ) 



W{-i'j'L'S' -^ -fJLS) 
4 (wa)3 



3 2J' 

fij'J'L'S' ^ 
2lo 1 
3" 2 J' + 1 



-\{-lJLS\\d\WJ'L'S')[^ 

^ jJLS) = 
\{jJLS\\d\WJ'L'S')\', 



(9) 



(10) 



where 7 denotes all quantum numbers other than J, L, 
and S. uj = E-^i ji — Ejj is the transition frequency from 
the initial state to the final state. With this definition the 
oscillator strength is positive for absorption and negative 
for emission. 

Using Eqs. ^ and pO|) and knowing El transition 
amplitudes between different states, we were able to cal- 
culate rates and oscillator strengths for the transitions 
involving the 5s^ ^Sq and the 5s5p ^Pq states. In Ta- 
ble |IV] we list the transition rates and oscillator strengths 
for the strongest transitions from the mentioned states. 
Note that the transition rates Wn-,iSo,^P° ^^'^ the os- 
cillator strengths /iSo,3po^„ were calculated with use of 
the theoretical energy levels. Where available we com- 
pare our results with other experimental and theoretical 



TABLE IV: Transition rates W„_>o and oscillator strengths fo^n for relevant energy levels in Sr. The results are compared 
with other available experimental and theoretical data. 







W^„^o(xloS-i) 




fo^n 


(0| 


\n) 


This work Other data 


This work 


Other data " 


5s' 'So 


5s5p 'P° 


186.0" 190.01(14) " 
191.6(1.1) -^ 
215 " 


1.82" 


1.92(6) 




5s6p iPi° 


1.49 1.87(26) " 
3.79 ' 


0.0058 


0.0072(10) 




5s7p 'Pi° 


5.13 5.32(61) " 
3.19 " 


0.15 


0.16(2) 


5s5p 'PS 


5s4d ■^'Di 


0.29 


0.088 






5s6s ^Si 


8.39 7.32 ■= 


0.173 






5s5d ^Di 


38.1 30.7 ■= 


0.395 






5p2 «Pi 


41.3 


0.418 






5s7s ^51 


2.28 1.80 " 


0.019 






5s6d ^Di 


14.3 


0.099 





"Parkinson ei al. 43](exp.). 

'These values are presented only for comparison. In all calculations performed in this work involving (^5o||d||^Pi°), we used 

the value of this ME obtained from Ref. [la |. 

■^Yasuda et aZ. JlpKexp.). 

''Nagel et al. [Ijfexp.). 

''These numbers were obtained from the values given in Werij et al. [4l[ with use of Eqs. ((9)1 and (|ll|l (see Subsection IIII Dl for 

details). 



values. As is seen from Table IIVI there is a reasonable 
agreement between the results of this work and other 
data. 

In certain cases we used for comparison the non- 
relativistic values of transition rates given in the litera- 
ture. In the LS coupling approximation, there is a simple 
relation between relativistic and non-relativistic reduced 
matrix elements (MEs) of the operator d. Since the op- 
erator d commutes with S we obtain [401, 



relation [4C 



{^JLS\\d\\-f'J'L'S') = 5ssV(2J+l)(2J' + l) 

s+L+j'+i I L J S 
J' V 1 



X (-1) 



{^LS\\d\WL'S).{ll) 



Knowing a non-relativistic transition rate we were able 
to determine the corresponding ME of the electric dipole 
operator {'^LS\\d\\^'L' S). If the LS coupling approxima- 
tion is valid, using Eq. ([TT|) . the non-relativistic ME can 
be related to the relativistic ME {-fJLS\\d\\i J' L' S'). Wc 
also need to account for the fact that in a non-relativistic 
case the transition frequency ujls between two states is 
given b y expression u^s = Eyj'L'S' - E^jls, where 
E-fjLS is the center of gravity of the respective multi- 
plet. For this reason, in general, wls can slightly differ 
from Lo — E^iji — E^j. Finally, using Eq. ^ we can find 
the relativistic transition rate. We used this approach to 
compare the relativistic transition rates obtained in this 
work with the non-relativistic values presented in [4l| . 

If the LS coupling breaks down, Eq. pT|) is no longer 
valid. In this case, to compare the relativistic transition 
rates given by Eq. ^ with the non-relativistic transition 
rates W{'j'L'S' -^ jLS), one should use a more general 



W{-/'L'S' -^ -/LS) 



(2L' 



"^{2J' + l)W{j'J'L'S' 



+ 1){2S'- 
jJLS). 



(12) 



JJ' 



In the right-hand side of this equation the summation 
goes over all possible values of J' and J. Consequently, 
we need to find all transitions rates (permitted by selec- 
tion rules) from the fine structure levels of one multiplet 
to the fine structure levels of another multiplet. 

To provide a straightforward comparison of the calcu- 
lation here with experimental data, lifetimes of the four 
states which dominate the ^P^ polarizability contribu- 
tions have been evaluated. For these four states, decay to 
the 5s5p ^Pj states is the only significant radiative decay 
channel so the lifetimes can provide direct information 
on the relevant matrix elements. A number of lifetime 
and transition rate measurements are available for com- 
parison [HI, m, m, m, 113, H H IM iH, in , however in 
many instances with limited accuracy. Table |V] summa- 
rizes the results where the 5s4d ^Di, 5s6s ^Si, 5s5d '^Di, 
and 5p^ "^Pi lifetime calculations are compared to avail- 
able measurements. In some cases, the measured life- 
times were reported for a particular J value in the ex- 
cited state multiplet, and in others only a mean lifetime 
for the entire multiplet was given, leading to complica- 
tions in the analysis. 

In the case of the 5s4d ^D state we also evaluated the 
total transition rate for the multiplet since to our knowl- 
edge the only lifetime measurements for the 5s4d ^D lev- 
els were performed on the entire multiplet. In the frame- 



work described above we found the W{5s4d^Dji -^ 
5s5p ^Pj) transition rates for all possible J and J', and 
using Eq. (fT2|) . we obtained, 



^y{2J'+l)W{bsAd^D.,. 
15 ^ — ' 
JJ' 



5s5p 3PJ) 



0.41xl0^s" 



in agreement with the experimental value W{hs4:d ^D -^ 
5s5p 3p°) = 0.345(24) x 10^ s'^ ^. 

The lifetimes of the 5s5d ^Dj states have been mea- 
sured with high accuracy (0.1%) in Ref. [50] and are in 
good agreement with other measurements as well as our 
calculated values. Ref. [50| also reported a 1% measure- 
ment of the 5p^ ^P2 lifetime, while other direct measure- 
ments of this multiplet yield consistent ^Pi lifetimes hav- 
ing accuracies at the 15-20 % level. Relative to 5s5(i ^Dj 
and 5p^ ^P.J: the 5s6s ■^S'l experimental data has larger 
scatters between different measurements. Notably all of 
these measured lifetime values agree well with our calcu- 
lations. 

Given the results in Table |V] and Table Hill the 
5s4d ^Di state should have the highest measurement pri- 
ority as it dominates the '^Pg static polarizability. It also 
has a large disagreement between the experiments |42l | 
and 545|, meriting further experimental investigations. 
The next priority goes to the 5s6s ^Si state due to the 
scatter in existing data and its large contribution to the 
a.c. polarizability of ^P§- Perhaps a good strategy is to 
measure its decay directly to individual 5s5p '^PJ states. 
Confirmation of the high accuracy result of Ref. [53| for 
the 5s5(i '^Di state and improvement upon the 5p^ ^Pi 
result listed in Table |Vl or alternatively, the use of the 
measured magic wavelength and clock laser light shift, 
can then be sufficient to determine the ^Pq polarizability 
at the 0.1% level. 



IV. CONCLUSION 

In this work we have carried out detailed calculations 
in response to the goal of further improving the accuracy 
of the Sr atomic clock to 1x10""'^^ and below. To focus 
on the outstanding problem of BBR-related frequency 
shifts, we calculated a.c. polarizabilities of the ^^o and 
■^Pg clock states. We verify our calculations with avail- 
able experimental data. For example, the theoretically 
calculated magic wavelengths for the ^S'g — "^Pg and the 
^S'g — ^Pi transitions are in 1% agreement with experi- 
ments. The agreement between theory and experiment 
on the a.c. Stark shift of the clock transition itself is also 
good. We have calculated individual contributions of six 
lowest-lying even-parity states to the polarizability of the 
■^Pg state aX oj — and for the wavelengths A — 698.4 nm 
and 813.4 nm. We determined four even-parity states 
whose total contribution to the static polarizability of 
the ■^Pg clock state is ~ 90%. Using the modern meth- 
ods of atomic calculations we can find the contribution 
of all the other discrete and continuum states (constitut- 
ing 10%) to the 3Po° polarizability at the level of 1-2%. 



TABLE V; Experimental lifetime data for the low-lying even- 
parity states. Theory values calculated in this work are pro- 
vided for comparison. ^D denotes measurement of the "^D 
manifold without definitive angular momenta. 



Excited state 


T (ns) 


This work 


5s4d ^Di 




2040 


3D 


2900(200) " 
4100(600) " 


2400 *" 


5s6s =^51 


15.0(8) "^ 

10.9(1.1) " 

12.9(7) ^ 


14.1 


5s5d ^Di 


16.49(10) 


16.9 


'D2 


16.34(13) » 


16.8 


'D, 


16.29(24) » 




w 


17.1(8) '^ 

16.0(6) '^ 

16.7(1.0) ■' 




5p^ ^Po 




7.8 


'Pi 


8.3(0.4) ■* 
8.8(1.2) ' 
10.2(2.4) ' 


7.6 


'P2 


7.89(05) » 
7.8(1.8) ■= 
8.3(0.4) "^ 


7.9 



"Reference | 42j . 

''This value is calculated using Eq. I|12| l. 

"^ Reference 

''Reference 

■^Reference 

■'Reference 

^Reference 

''Reference 

'Reference 



For this reason, if the contributions of the four states 
identified here are experimentally determined with 0.1% 
accuracy, the same level of accuracy can be obtained for 
the total polarizability of ^Pq. In the near future we 
plan to undertake experimental measurements to deter- 
mine the oscillator strengths for the four identified states. 
Measurements could include transition linewidths, power 
broadening coefficients, direct lifetime determinations of 
individual J levels, and improved determination of the 
clock laser a.c. Stark shifts. These experimental mea- 
surements can be further combined with the well-known 
value of Ag. The experimentally determined values will 
be used to refine the theory calculations presented here 
to reach the goal of determining the polarizability of ^Pq 
at 0.1%. 
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